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Abstract 

Some properties of small and medium sodium clusters are described within the RPA approach 
using a projected spherical single particle basis. The oscillator strengths calculated with a Schiff- 
like dipole transition operator and folded with Lorentzian functions are used to calculate the 
photoabsorbtion cross section spectra. The results are further employed to establish the dependence 
of the plasmon frequency on the number of cluster components. Static electric polarizabilities of 
the clusters excited in a RPA dipole state are also calculated. Comparison of our results with the 
corresponding experimental data show an overall good agreement. 

PACS numbers: 36.40.Gk, 36.40.Vz, 36.40.Ei 
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I. INTRODUCTION 



Although the field of metallic clusters is very old, it became very attractive for both 
theoreticians and experimentalists only since 1984, when the pioneering paper of Knight et 
al [l| pointed out the electronic shells in alkali-metallic clusters. Some notable contributions 
in the new era of the field have been reviewed by several authors 

aaa. 

The clusters which are most interesting to be studied seem to be those^with a moderate 
number of atoms. Indeed, 



quantum-chemical methods 



'or such systems neither statistical models [s] nor ah — initio 
6j are justified. Instead, the mean field approach is vastly used. 
Several solutions defining the mean field for the single particle motion have been employed 
along the years. Among them, three procedures are to be distinguished: i) solving the 
Kohn-Sham equations {7], ii) assuming that the positive charge of the ionic core is uniformly 
distributed in a sphere of radius R. This is known in the literature as jellium hypothesis 
js, 9]. iii) postulating the average potential [lo| . 

The researchers in this field payed attention especially to spherical clusters since to this 
shape are associated notions like shell structure or magic numbers. However, there are 



some features like the detai 
plasmon energies Ijl, Q 



ed structure of the abundance spectrum jll], or the split of 



14j, which cannot be explained assuming a spherical symmetry 



for the mean field. The first paper devoted to the deformed clusters was due to Clemenger 
and published in 1985 [la ]. The author adapts the Nilsson model formulated for nuclear 
systems 16], by ignoring the spin-orbit term. The resulting model is, therefore, referred to 
as Clemenger-Nilsson (CN) model. The CN model was very successful in explaining several 
properties which depend on the cluster shape and which could not be described within a 
formalism using a mean field with spherical symmetry. The model is suitable to describe 
the single particle properties in the intrinsic frame and especially for strong coupling regime 
when the wave function of the whole system can be factorized into an intrinsic part and a 
Wigner function accounting for the rotational degrees of freedom. For clusters exhibiting an 
axial symmetry the projection of angular momentum on the symmetry axis, denoted by K 
is a good quantum number. Of course the angular momentum itself is not a good quantum 
number. Having in view the fact that the measurements are achieved in the laboratory 
frame where the rotation symmetry is valid we have to pay attention to this feature. Indeed, 
there are many properties which are very sensitive to the change of angular momentum of 
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the system. Moreover in most cases "K" is not a good quantum number and therefore the 
factorization mentioned above is not possible. The typical case of this kind is that of systems 
of triaxial shape. A many body treatment of such situations would require a subsequent 
projection of the angular momentum. Such an operation is technically very difficult to be 
achieved and to our knowledge up to now only approximate solutions have been adopted. 
In a previous publication, one of us (A. A. R.) proposed a solution for construc ting a single 
particle basis with good spherical symmetry and depending on deformation [17|. In the 
quoted work the model ability to account in a realistic fashion for the main features of the 
deformed clusters has been successfully tested. Thus, the cluster shape, the magic numbers, 
and super-shell effects have been determined and a good agreement with the data as well as 
with the previous theoretical results have been obtained. 

In a subsequent paper we continued the exploration of the deformed cluster properties, 
within the projected spherical basis introduced in ref. The single particle symmetry is 
actually the mean field symmetry. This feature is an important property related with the 
charge distribution of the valence electrons. Such a structure becomes very important when 
one studies the response of the cluster to the action of an external electromagnetic field. 



Guided by a possible paralle 
out in Ref. 
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ism between atomic clusters and nuclear systems pointed 
2l| . in Ref. |22j one of us (A.A.R.) studied, in collaboration, properties 
like skin structure, empty center, hard center, cluster subsystems and halo behavior. Such 
properties have been seen, indeed, in the structure of the charge density function. Based on 
phenomenological arguments, quantitative results for polarizability and plasmon frequency 
have been derived. 

In this paper we study the collective many body properties of light and medium sodium 
clusters. The paper is structured as follows. In Section II we briefly review the projected 
single particle basis. In Section III the basic equations specific to the particle-hole (ph) 
Random Phase Approximation (RPA) are written down. The RPA wave functions are used 
in Section IV, to treat the electric dipole transitions and the photoabsorbtion cross section. 
In the framework of phRPA the expressions for system polarizability are analytically derived, 
in Section V. Numerical applications for Na clusters with the number of components varying 
from six to forty have been performed and the results are described in Section VI. The final 
conclusions are summarized in Section VII. 
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II. PROJECTED SPHERICAL SINGLE PARTICLE BASIS 



We restrict our considerations to the energy domain of laser beam experiments (i.e. opti- 
cal domain) where only the valence electrons may be excited and de-localized, the remaining 
ones defining the atomic core. Under these circumstances we could study those atomic clus- 
ter properties which are mainly determined by the valence electrons. The picture is even 
more simplified if the cluster building block is an alkali-metal. In this case the interacting 
system of electrons and positively charged ions is replaced by a system of interacting elec- 
trons moving in a mean field which accounts for the influence of the ionic core on the single 
particle motion. Each atom of a given cluster is represented by one valence electron. 

In a previous publication the mean field for the valence electrons was defined with 
the help of a model Hamiltonian associated to the particle-core interacting system. 



A=0,2 -A<^t<A 

= H p + H c + H pc , (2.1) 

where are A-pole shape variables defining the deformed ionic core through the surface 
equation: 



R = R (l+J2 E a V> 

V A=0,2 -A<^t<A 




(2.2) 



The volume conservation condition allows us to relate the monopole and quadrupole coor- 
dinates: 



1= E i^i 2 - ( 2 - 3 ) 



'Air 

-2<A*<2 



For what follows it is convenient to introduce the boson operators b\ 



(bl + (-Y%-») ; tt 2 , = 1= (-b 2 » + • (2-4) 



ky/2 

The core subsystem is considered to be a harmonic quadrupole boson Hamiltonian: 



H c = cu c J2bt^. (2.5) 
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Let us consider the coherent state 

\V c ) = exp[d(bt -b 20 }] |0>, (2.6) 

where |0) denotes the vacuum state of quadrupole bosons. This wave function is a coherent 
state: 

b 2m \^c) = <W|*c) (2.7) 

and describes the ground state of a deformed quadrupole boson Hamiltonian. Moreover, the 
expected value of the static quadrupole moment in the state |^/ c ), is proportional to d. Due 
to this feature, d is referred to as the deformation parameter. 

Note that averaging H with |\I/ C ), we obtain a deformed single particle mean field which 



is similar to that used by Clemenger in Ref. [15|. On the other hand, averaging H with 
an eigenstate \nlm) of H p , one arrives at a deformed quadrupole boson Hamiltonian which 
admits |\l/ c ) as ground state if a suitable deformation parameter is chosen. These properties 
suggest that the particle-core system wave function might be described by a transformation 
applied to the product function with the factors \nlm) and |^/ c ). 

The deformation of the single particle mean field in the Clemenger's model, 5, might be 
related to the deformation parameter d involved in the coherent state defined above. Aiming 
at this goal, we require that the average of H pc with |\l/ c ) is identical to the deformed single- 
particle potential from the Clemenger's Hamiltonian. This supplies us with the relation 

| = ^(«i-°3, («) 

where fl z and Q±_ denote the frequencies along and perpendicular to the symmetry axis, 
respectively. This relation yields a simple equation for the two deformation parameters: 

0.693fc5 = d. (2.9) 

The mean field defines an intrinsic frame of reference for electrons. In the laboratory 
frame the system is described by states having good angular momentum, due to the rotation 
symmetry of the model Hamiltonian. A basis in the laboratory frame for the composite 
system, of electrons and core, may be obtained by diagonalizing H in a particle-core product 
basis, with components of definite angular momenta. However, to use such a basis in a 
Random Phase Approximation (RPA) formalism is quite a tedious task. 
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A great simplification is obtained if, instead, we use a projected spherical single particle 
basis. Obviously, even if the deformed set of generating functions is orthogonal the angular 
momentum projected set is not orthogonal. Fortunately, an orthogonal set of projected 
states is obtained if we chose an appropriate single particle state factor. 

Thus, it can be proved that the following subset of projected states is orthogonal: 



<I>im;v (nl; d) = Aft(d) [Pl H \nlI)y c {d)} X *, for I ^ 0, I 
0oo ; <x(nZ;d) = N^{d) [p ° [\nl)s] l+lfl ^ c (d)} X o 



even, 



for / = 0, / = odd, 



(2.10) 



where s denotes the spin operator and Xa is the bi-spinor component. The standard notation 
for the angular momentum projection operator was used: 

21 + 1 



P 



MK 



D\* K (n)R(n)dn. 



(2.11) 



The norms Ml^d) of these projected states are 



ui\ 2 ( AT-(C) 

j 



, for / ^ 0, I — even, 
for / = 0, I = odd, 



(2.12) 



with Nj^ denoting the norm of the J component projeted from the deformed state \I/ C ), 
describing the core. If one neglects the matrix elements with Al = ±2 and An = 2, the 
eigenvalues of H within the projected spherical basis can be fairly well approximated by the 
average values: 



e^d) = < 4>iM;a( nl 'i d)\H\(f) IM . a (nl; d) >= Huj ( N + 



D 



1(1 + 1)- 



N(N + 3) 



+ 



' L v 



(2.13) 



The expected va 
analytically in Ref. 



ue of the boson number operator as well as the factor Fn were calculated 



~v\ . The energies e^(c?) depend on d in the same manner as the single 
particle energies in the Clemenger's model depend on the deformation 5. This comparison 
has been performed in Ref. 17|. The deformation parameter is fixed so that the cluster 
energy is minimum against any variation of d. 
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As shown in Ref. [3], although the functions are characterizing a particle-core system, 
they can be used to describe a many fermion system. Indeed, in calculating the matrix 
elements of an one body operator, one integrates first on the core coordinates and then on 
single particle coordinates. Finally, the matrix element is written in a factorized form, one 
factor describing the matrix elements between spherical wave functions and the other one 
carrying the dependence on the deformation. Moreover, due to the specific properties of the 
coherent state the projected spherical states can be used for calculating the matrix elements 
of a two body operator. Indeed, in Ref.[19| we have proved that the matrix elements of a 
two body interaction are practically equal to the matrix elements between states projected 
from the product of two particle states and a common core's coherent state. 

To conclude the single particle energies are approximated by those given in Eq. (2.13) 
while the corresponding wave functions by Eqs. (2.10), (2.12). Actually, these are input 
data for the treatment of a many body Hamiltonian associated to a system of interacting 
valence electrons moving in the mean field presented above. 



III. RPA DESCRIPTION OF THE COLLECTIVE DIPOLE STATES 

We assume that the valence electrons moving in the mean field of the ionic core and inter- 
acting among themselves through a Coulomb force are described by the model Hamiltonian 

H = H mf + V c , (3.1) 

which is a sum of the static self-consistent, single-particle Hamiltonian H m f and of the two- 
body residual Coulomb interaction Vc- The residual two-body interaction is given, in the 
local-density approximation, by 

V(n - f 2 ) = + d ^f^5{r 1 - f 2 ). (3.2) 

\r 1 -r 2 \ dp 

The quantity Vxc = dexc[p]/dp is the exchange-correlation potential in the ground state. 



We use the exchange-correlation energy density Exc of Gunnarsson and Lundqvist 
in Refs. 



24 



23| as 



251 ]. Thus, the following expression for the exchange-correlation potential Vxc 



given in atomic units, is obtained 



1 222 / 11 4 \ 

Vxc(r) = -- -- - 0.0666 In 1 + — - . (3.3) 
r s [f) \ r s (r) J 
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Here r s (f) = [3/47rp(r)] 1//3 is the local value of the Wigner-Seitz radius. The two body 
interaction is expanded in multipoles, the A-pole term having the expression {24]: 

1 



V(r h r 2 ; A) 



.2 r < 



A + dV xc [p]5(r 1 -r 2 )2\ 



„A+1 

> 



dp 



Art 



(3.4) 



where r< = minfri, r 2 ) and r> = maxfri, r 2 ) 



For A = 1, the two body interaction has the expression: 



F(p,r s ) 



e 2 ^ + F(p,r s ) 



5{r x - r 2 ) 



4:iipr s 



1.222 + 



0.759r, 



r s + 11.4 



(3.5) 



The one body term represented by the mean field H m f and the A = 1-pole term of the 
two body interaction are treated within the RPA formalism which defines an operator 



C U = [ X ph(cl c h)^ ~ Y ph(c[c p ) llM 

ph 



(3.6) 



subject to the restrictions: 
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Written in a matricial form, the RPA equations provided by 13.71 look like: 

A B ){ x )-4 X " 

-A* -b* My I \ Y n 

The sub-matrices A and B have the expressions: 



A(ph;p'h' 



(e p - Sh)8 p ,pi8 hth > - B(ph;p'h'), 
tip 

y 



Birlr.,/1/) 2 IpIhlhI :' Ih ' lh ' ^^oo^Vfc^^^^^^fi^jV) 



where 



(3.7) 
(3.8) 



(3.9) 



(3.10) 



R(pti; hp') = J rldr l rldr 2 n np i p {ri)n nhl i h Xr2)V{r u r2] l)7^ fe (n)^ V J p ,( r2 )> (3.11) 

with IZn^ir) being the radial part of the single-particle projected wave functions. Since 
the two-body interaction consists of two terms, the Coulomb and the exchange term, corre- 
spondingly the factor R(ph'; hp') split into two parts which are given in Appendix A. The 
factor 2 of the sub-matrix B accounts for the spin degeneracy. 



8 



The factor f(d) carrying the dependence on deformation parameter d, involved in 
Eq. fl3.10l) . have the expressions given in Appendix A. Eq. (l3.9p determine the amplitudes 
X and Y up to a multiplicative constant which is fixed by the normalization condition (13.81) : 

£[W- W] =i- ( 3 - 12 ) 

ph 

The compatibility condition for the set of equations (I3.9P is an equation for u. The number 
of solutions, for the RPA equations, is equal to the number of dipole particle-hole (ph) 
configurations, hereafter denoted by N s . To the solution u k , the amplitudes X^ h and Y* h 
correspond. These amplitudes characterize the phonon operator ClJk) which may excite 
the cluster ground state |0) to a one phonon state: 



CV(fc)|0) = 0, \l k p) = Cl(k)\0). (3.13) 
There are several procedures to solve numerically the RPA equations. Here we have 



26|. 



adopted the method proposed by Rowe in Ref. 

Concerning the RPA description, we would like to comment on the following features: 

a) Due to its specific structure, the state with 1 = and / = odd may be related, by 
the ph dipole operator, to states with either the spin up or the spin down. Consequently 
the RPA matrix has not a block structure, each block being characterized by an unique 
orientation of the electron spin. The coupling terms are however small and bring negligible 
contribution. For this reason we ignore from the beginning the terms which flip the spin of 
the ph matrix elements. 

b) In the definition of the phonon operator (13. 6p the summation involves both states (2.10) 
with the factor state X1/2 and with the bi-spinor X-i/2- We, conventionally, call the resulting 
operator as the extended phonon operator. In virtue of a) we may restrict the summation 
to one component and therefore work with the so called reduced phonon operator. The 
normalization to unity of the two phonon operators suggest a simple relationship between 
their defining amplitudes. 

In our approach the reduced phonon operator has been used, otherwise the degeneracy 
of the states with spin up and with spin down has been carefully implemented whenever the 
matrix elements between RPA states were calculated. 
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IV. El TRANSITIONS AND PHOTO ABSORBTION CROSS SECTION 



The reduced probability for the dipole transition |0) — > |1+) can be written as W. 

b(ei,o+^i+) = \(o\\M(ei)\\i+)\ 2 , (4.i) 

where 

(0||A<(JBl)||l+> =J2Up\\M(El)\\h) [X; h + (-y* +I %l] (4.2) 

ph 

are the reduced matrix elements of the dipole operator A4(E1), between the specified RPA 
state. The one phonon state |l n ) is characterized by the RPA amplitudes X™ h and Yp h , 
obtained by solving Eq.(3.7). 

Instead of the usual transition dipole operator, we make use of a modified operator which 



is similar to the so-called Schiff moment from nuclear physics 



3l| 



Mm = ^ Y ey 1 ,(Q) [r - --J . (4.3) 

Here r s is the Wigner-Seitz radius and is equal to 3.93 (a.u.) for Na clusters. The corrective 
component, involved in the dipole operator, relates particle and hole states characterized 
by AN = 3, which results in modifying the strength distribution among the RPA states. 
Such an effect is obtained in a natural manner, i.e. using the standard form for the dipole 
transition operator, if the mean field potential for the single particle motion involves higher 
powers of the radial coordinate. 

We define the oscillator strength f n per atom as 

In , 

where u n is the RPA excitation energy corresponding to the n-th order solution of RPA 
equation, and 

S(£l) = 5> n £?(M,0 + ->l+), (4.5) 

n 

such that J2 n fn = 1- 

To calculate the photoabsorbtion cross section per atom, <j(uj), one folds the oscillator 
strengths, which are just vertical straight lines, with Lorentzian shapes normalized to unity 



[1] Throughout this paper the Rose's convention for the reduced matrix elements are used 
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as follows: 

a(m) = C ^ / n £(^; Wn, T„) , (4.6) 

n 

where T n denotes the full widths at half maximum of the Lorentzian profiles, and is provided 
by fixing the damping factor 7 = V /u r with uj r being the energy of the resonance peak. This 
damping factor varies, in our calculation, in the range of 0.06-0.135 [3| which is appropriate 
for the room temperature. Indeed, the damping factor is considered to be caused by the 
coupling of the electronic dipole oscillation to the thermal fluctuations of the cluster surface 



32| | . The proportionality coefficient C is given by 



C=^^ = 1.0975 (eVA 2 ). (4.7) 
m e c 

This value is obtained by normalizing the photoabsorbtion cross section such that the area 



per de-localized electron under the photoabsorbtion curve is constant 



35|, 



DC 







a(u)du = , 4.8 

m P c 



which is consistent with the value of the dipole sum rule. Note that the results for the 
reduced El transition probability as well as for the photoabsorbtion cross section depend 
on the single particle features specified by the mean field parameter Huq. The oscillator 
parameter of energy quanta hujQ, depends on the cluster atoms number as Ep/ftf 1 / 3 , where 
Ep is the Fermi energy, which is about 3 eV for a spherical cluster. However, in general, Ep 
has not a constant value for all clusters. Keeping this in mind, we use a Fermi energy Ep(N) 
which also depends on the number of cluster's atoms. The Af dependence is extracted by 
interpolating the results provided by a least square fit for the experimental photoabsorbtion 
cross section spectrum. In our calculations one also needs to know the oscillator length 
b. This depends on the choice of oscillator energy quanta, and therefore exhibits the M 
dependence given by b= [% 2 /m e E F {U)\ ll2 N l ' & . 

V. ELECTRIC POL ARIZ ABILITY 

In the classical picture, the static electric polarizability of a jellium metal sphere of radius 
R has the expression: 

a = R 3 , (5.1) 
where R, given in terms of the Wigner-Seitz radius r s , is R = r s A^ 1//3 . 
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Quantum mechanical effects determine corrections to the classical results for plasmon 
energy and polarizabilities. The electron density is not going sharply to zero at the cluster 
surface but reduces gradually at the surface and moreover extends significantly beyond to 
jellium edge. The spill-out electrons produce a screening effect against external fields which 
results in changing the classical result for polarizability to: 

a = (R + A) 3 . (5.2) 

The radius shift A can be expressed in terms of the fraction of the total number of electrons 
which are spilled-out the jellium sphere and the final result for the static polarizability reads: 

« = R 3 (l + ^f) , (5-3) 

where Af sp denotes the number of spilled-out electrons. 

We propose a method for obtaining the number of spill-out electrons by means of the 
RPA formalism, using the RPA eigenstates for computing the average value of number of 
particles operator N = X^ c ! c «> where the summation run over all particle and hole states. 
Thus the operator to be averaged is given by 

# = £<&+£ c U- ( 5 - 4 ) 

p h 

The averaging operation will be constrained by the condition that the radial integrals from 
all scalar products involved will have the limits [R, oo), instead of [0, oo), where R = TgH 1 ^ 
is the radius of the metal sphere which is supposed to contain all electrons. In this manner 
we will get the average number of electrons which are beyond this sphere. The particle 
number operator can be written in a second quantization like form 

kk' 

Averaging N with the RPA eigenstates we get: 

(N) = E(^l^l^) = 2 EE [H^KJpY + WY p k XY] • (5.6) 

k k ph 

Here, the overlap of two single-particle projected functions corresponding to two particles 
or two holes states are expressed through the product of the statistical factor and the 
radial integral 

POO 

/* = / [n nih (r)] 2 r 2 dr, t = p,h. (5.7) 
Jr 
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VI. NUMERICAL RESULTS 



First, we identify the dependence of the Fermi energy on the number of atoms in cluster. 
The adopted procedure is as follows: For each cluster we determine the Fermi energy which 
corresponds to the best agreement of the calculated photoabsorbtion curve with the exper- 
imental points. Further, the obtained values are interpolated with a third order polynomial 
in A/" 1//3 (see Fig{T]). The polynomial obtained in this way determines a Fermi energy varying 
in the range of 3.3-3.75 eV. The RPA calculations make use of the Fermi energies lying on the 
interpolating curve. In this way the RPA results depend on the number of atoms, by means 
of the oscillator energy quanta hu = Ef{M)/M 1 ^ and the oscillator length b = ^h/m e uj . 
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FIG. 1: The Fermi energies (black circles) which reproduce the experimental photoabsorbtion cross 
section spectra, are interpolated with a third order polynomial in AA 1//3 : Ep{M~) = — 0.886466AA + 
7.47916AA 2 / 3 - 20.3871AA 1 / 3 + 21.4339. 

The RPA calculations are performed for the restricted subspace of the particle-hole exci- 
tations, including only the AN = 1 and AN = 3 transitions. We use the method described 



in Ref . 



261 ] to solve the RPA equations with real sub-matrices A(ph;p'h') and B(ph;p'h'). 



The single particle energies £j involved in the expression of the sub-matrix A(ph]p'h') are 
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provided by the formula (2.13) in units of huio. Once the RPA amplitudes and energies are 
determined, the electric dipole transition probabilities and the oscillator strengths can be 
computed by means of Eqs. (4.1) and (4.4). The biggest values for the transition prob- 
abilities define the collective states. The corresponding RPA energies are located mainly 
in two regions associated to the first collective dipole state, with the major contribution 
coming from the AN = 1 excitations, and the second collective dipole state, which is due 
to the AiV = 3 transitions. The energy of the state characterized by a dominant transition 
probability varies from one cluster to another. Thus, the energies of the first and second 
collective states, depend on A/" 2 / 3 and A/" 1 / 3 respectively, as shown in FigJ2 Therefore, the 
first collective state seems to have a surface mode behavior, while the second collective ex- 
citation exhibits both a surface and a volume feature. In a short interval of M the volume 
character of the collective dipole mode of larger energy, prevails. 



1 1 1 


1 1 1 1 1 




; -v 










— •— Energy of the first 




collective dipole state 


■ a) 

1 . 1 
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E =0.91 046N 2/3 -5. 4091 5N 1/3 +1 3. 9351 



FIG. 2: The RPA energy of the first collective dipole state decrease very slowly and almost linearly 
in AA 2 / 3 (a). For the second collective dipole state the energy is decreasing faster by means of 
a parabolic law in AA 1 / 3 (b). Theoretical results for all the clusters in the range of 8-40 atoms 
constituents, are presented. 

In Table I, we collected the RPA amplitudes X p h and Y p h associated to the first collective 
dipole states in some Na clusters. From there one sees that notable contributions are brought 
by several particle-hole configurations, which contrasts the situation of the non- collective 
states where a certain component is by far dominant, while the others are negligible. 

In order to compare the predictions of the present formalism with the experimental results 
for photoabsorbtion cross sections, we have also to take into account the coupling of the 
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electronic dipole oscillations with the thermal fluctuations of the cluster surface [321 . |34j . In 
the generalized picture, the photoabsorbtion cross section per atom for a spherical cluster 
that is much smaller than the photon wavelength [33J , have the plasmon energy dependence 

a(cu) = An- — 2 , (6.1) 



where T is an averaging parameter [32J. Here the surface plasmon pole approximation has 
been adopted. The resonance energy was denoted by uj t while V is just the resonance width. 
Keeping close to the general description of the photoabsorbtion cross section, the above 
formula can be simulated approximately by folding one oscillator strength line, predicted by 
the RPA calculation with a Lorentzian function displaying a damping factor 7 = Y/u r ~ 0.1, 
which is appropriate for room temperature. This procedure is extended to folding all the 
strengths provided by the RPA calculation, according to Eq.(4.6). Indeed, to each energy 
ojk one associates a Lorenzian centered in Uk and having a width given by the product of a 
constant damping factor 7 and the Lorenzian resonance energy. In Figj3] theoretical curves 
of photoabsorbtion cross section per atom are displayed versus the wavelength, for several 
clusters for which experimental data are known in a large interval. Note that the RPA 
energies were given in terms of the corresponding wavelength A = lnTicjiOk- Apart from 
some discrepancies in the form of the curves, like plateaus or number of peaks, the present 
RPA calculations provide an overall good agreement with the experimental measurements 
in the visible range of energies. An interesting shape of the photoabsorbtion cross section is 
obtained for Nan cluster, which has a triaxial shape reflected in the fragmentation of the 
oscillator strength spectrum into three comparable peaks. These three distinct peaks are 
identified as plasmon frequencies associated to those three axis which define the triaxiality of 
the cluster. In a spherical cluster like Nag and Na2o these three frequencies are degenerate 
which results in having only one peak in the photoabsorbtion cross section curves. Also, our 
calculation predicts the double peak structure of photoabsorbtion spectra for clusters Na w , 
Nan and iVa 12 . Actually this structure suggests an axially symmetric shape. 

Taking a look at the photoabsorbtion cross section spectra in the ultraviolet region of 
wavelengths (Fig. HI), where the second collective dipole state shows up, it is noticed that 
the corresponding photoabsorbtion peaks are also fragmented for some clusters. For some 
clusters there exists only one peak in this wavelength domain and this happens due to the 
near spherical form of these clusters. Indeed, the single peak shape is caused by the fact that 
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for spherical clusters there is only one collective state, while for strongly deformed clusters 

this state is split into two parts, one consisting in two degenerate states of energy ui x = lu v , 

PI 

and another one of energy uj z [22|] . The ordering of the two energies depends on the cluster 
shape. For prolate clusters uj z < u x , while for oblate clusters the ordering is changed. It 
is natural to suppose that the highest peak corresponds to the double degenerate energies 
uo x = uj y . Keeping this picture in mind and inspecting the crossections shown in FigJH 
one can conclude that for most clusters with the second collective dipole states fragmented, 
the ordering of the larger and the smaller peaks suggest a prolate shape. By contrast, the 
cluster iVai4 has an oblate shape. This reasoning agrees perfectly well with the experimental 
data. Also, it is known that Na$ is a spherical cluster, which is reflected in a single almost 
degenerate peak of the second collective dipole state. 

In Figj5] theoretical estimations for the static electric polarizability per atom, normalized 
to the value of polarizability of neutral Na atom are shown together with experimental 
values and LDA (Local Density Approximation) results, for a few clusters. The three data 
sets are compared with each other as well as with the bulk limit result, which is associated 
to the classical polarizability of a metal-sphere. An excellent agreement is obtained for 
most clusters' polarizabilities. However, there are few cases where big discrepancies are 
recorded, namely Nag, Naxo, Nan and Na^. The noteworthy fact is that these clusters are 
most fragmented, some of them displaying a smooth two peak structure or even continuous 
plateaus of the photoabsorbtion cross section spectrum. 
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FIG. 3: Experimental (circles) and calculated (solid lines) photoabsorbtion cross sections vs wave- 
ength for Na clusters of small and medium sizes. The experimental points are taken from ref. 
2?j], where the open circles corresponds to data taken with a flash-lamp laser and have about 20% 
statistical errors. Solid circles have smaller statistical errors, about 5 — 10%, for which cw lasers 
were used. The calculated curves are resulting from the folding of the RPA oscillator strengths 
with Lorentzian shapes normalized to unity. The Lorentzian shapes are specified by the relative 
width parameter 7 = T/tu r , which varies from tie cluster to another in the range of 0.06 to 0.135, 
uj r being the peak energy of a given Lorentzian. 
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FIG. 5: The predicted static electric polarizabilities per atom for Na_\f clusters normalized to the 



measured polarizability of Na atom (Ref. 29() (open 
experimental data (black triangles) from reference 
LDA approach. The bulk limit is also visualized. 



circles), are compared with the corresponding 



281 ] and with those given in Ref. 



30l | with a 



18 



TABLE I: The dominant RPA A^-amplitudes of the first collective dipole states achieved by 
AN = 1 transitions for some of the Na clusters. The corresponding dipole ph configuration as well 
the RPA energies are also mentioned. Note that the heavier are the clusters the more collective is 
the depicted RPA mode. Also due to the repulsive character of the two body interaction the order 
of the RPA root is increasing with M despite the fact that the energy is slightly decreasing. 
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VII. CONCLUSIONS 



The main results described in the previous Sections can be summarized as follows. The 
atomic clusters were replaced by a set of valence electrons moving in a mean field and in- 
teracting among themselves through a Coulomb and an exchange interaction. Although 
not presented analytically, the mean field is defined by a set of orthogonal projected spher- 
ical single particle states and a set of corresponding energies. The two-body interaction 
is expanded in multipoles, the expansion being truncated at A = 1. The mean field and 
the two-body dipole interaction are treated within the RPA approach which defines a set 
of particle-hole like phonon states. The El transition from the ground state to excited one 
phonon states were calculated. The collective El transitions to the low energy states, around 
2.5 eV, are determined by a coherent contribution of the AN = 1, ph configuration. There 
are also states, of energies about 5-6 eV, which are collectively populated with AN = 3 ph 
transitions. As shown in Fig. 2, the first collective mode has a surface character while the 
second one seems to be mainly of a volume type. 

The photoabsorbtion cross section was obtained by folding the El strengths carried by 
the RPA one phonon states, with Lorentzian with a damping factor varying from one cluster 
to another, in the range of 0.06-0.135. The result shown in Eq.(4.6) is plotted in Fig. 3 as 
a function of the wavelength associated to the RPA states, and compared with the experi- 
mental data. The comparison reveals a reasonable good agreement. For the light clusters, 
N=9,10,ll,12, the figures exhibit two peaks, which suggest the existence of two modes cor- 
responding to oscillations along and perpendicular to the symmetry axis, respectively. For 
N=14, the results of our calculation indicate a triaxial shape. Similar fragmentation is also 
obtained for the volume plasmon resonance, shown in Fig. 4. According to the relative mag- 
nitude of the two peaks, we concluded that the clusters with N=9,10,ll,12 has a deformed 
prolate shape while Nan seems to be of an oblate shape. Perhaps a smaller damping factor 
would provide evidence for a triaxial shape also for the volume mode in Nan. 

It is noteworthy the fact that the fragmentation effects seen in both low (surface like) and 
high energy (volume like) modes cannot be described with a spherical single particle basis. 

Based on calculations for the number of spilled-out electrons, we calculated the electric 
polarizability which is compared with the experimental data, the bulk limit result as well 
as with the LDA (local density approximation) predictions. Except for the clusters with 
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N=9,10,ll,12,24, the calculated polarizabilities agree quite well with the corresponding ex- 
perimental data. It is interesting to remark that concerning the clusters for which large 
discrepancies were recorded, some of predictions (N=9, 10,24) lie close to the LDA results 
while others (N=ll,12) agree with the bulk limit of Na clusters. 

The final conclusion is that the projected spherical single particle basis seems to be an 
useful tool for describing the many-body features of deformed atomic clusters. 
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VIII. APPENDIX A 

The two terms of the two body potential l / (r 1 ,r 2 ; 1) determine two following two body 
matrix elements: 



R(pti, hp') 
Ri(ph',hp') 



R 1 (pti,hp') + R 2 (pti,hp'), 



/ / r fRn P i P (r i )R nh i h (ri)dri R n > p i p (r 2 )R n > h i> h (r 2 )dr 2 , 



Jo Uo 





(A.l) 



The factors / involved in the RPA matrices have the expressions: 



x l p C\ h h lp +1 1; W{1 l p I h l h ; l p + 1 1), 




1 



(A.2) 



for (I p , lp) 



(Q,odd) and (Ihjh) ^ (0,odd), 




V2k + 3 



1 



(A.3) 



for (4, l h ) 



(0,odd) and (I p ,l p ) ^ (0,odd), 



Ktid) = v(lMh)K P p i P (dW% lh (d) 




(A.4) 
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for (Ip,l p ) 7^ (0, odd) and (IhJh) 7^ (0, odd), where v{Ii) = 2 j.+i° ^ s the statistical factor 
connected to the fact that the energy level with a given / have the degeneracy 21 + 1 and 
it contains also the I, —I degeneracy which add a factor 2 if J ^ 0. The Clebsch- Gordon 
coefficients appearing in the above expressions take proper care of the angular momentum 
coupling, as well as of the parity-conservation conditions. 
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